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Abstract 

We examine dynamic structure factors of spin-1/2 chains with nearest-neighbor interactions 
of XX and Dzyaloshinskii-Moriya type, and with periodic and random changes in the sign of 
these interactions. This special kind of inhomogeneity can be eliminated from the Hamiltonian by 
suitable transformation of the spin variables. As a result, the dynamic structure factors of periodic 
or random chains can be computed from those of the uniform chains. Using the exact analytical 
and precise numerical results available for the uniform systems we illustrate the effects of regular 
alternation or random disorder on dynamic structure factors of quantum spin chains. 
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I. INTRODUCTION. JORDAN- WIGNER FERMIONS AND DYNAMIC QUAN- 
TITIES 



Quantum spin chains have received much attention during the last more than 70 years 
for several reasons. First, they provide an excellent ground for studying rigorously quantum 
many-particle phenomena. Second, owing to the tremendous progress in material sciences 
(as well as the recent availability of optical lattices for trapping atoms in artificial crys- 
tals) many real-life systems, which can be modeled as quantum spin chains invented by 
theoreticians, have become available. That opens the possibility to compare the results of 
accurate theoretical calculations with experimental data. Dynamic quantities for quantum 
spin chains are of special interest and importance. On the one hand, their study, as a 
rule, is a harder problem in comparison with similar studies of static quantities. On the 
other hand, dynamic quantities are related to experimental data obtained in scattering and 
resonance experiments which yield valuable information about the magnetic structure of 
materials provided one has a reliable theory for their interpretation. Therefore, the theo- 
retical analysis of the dynamic quantities for quantum spin chains is significant both from 
theoretical/academic and experimental/practical points of view. 

Since the early 1930s the Bethe ansatz has been known as a powerful method of exploring 
quantum spin chains. However, only recently it has become possible to calculate quantities 
such as norms of and matrix elements between Bethe ansatz states which are necessary to 
calculate dynamic quantities. For recent Bethe ansatz results on the ground-state dynamic 



tne _grou 



structure factors of the spin- 1/2 XX Z Heisenberg chain see Refs. hi la la 

We also mention 

here the field-theoretical approaches for evaluation of the dynamic quantities 4( which, 
however, are restricted to low-energy physics only and therefore can only give the dynamic 
quantities in a small part of the plane wave-vector k - frequency u (hereinafter the k- 
uj plane). Traditionally, those calculations were performed for the Tomonaga-Luttinger 
model, which describes one- dimensional spinless fermions moving in a continuum, with linear 
dispersion relation. Recently, however, the curvature of the dispersion relation has been 
taken into account in calculating the properties of quantum wires js|. The spin-1/2 XX Z 
chain is a lattice system closely related to these continuum models. The low-energy and 
long-wavelength limit of its ground-state zz dynamic structure factor was recently studied 
by combining several analytic and numeric techniques Recently j/J an extension of the 
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density-matrix renormalization group (DMRG) method was proposed which allows for the 
calculation of real-time correlation functions of XX Z chains at arbitrary finite temperatures 
in the thermodynamic limit. However, numerical limitations presently restrict the time 
range over which results are reliable to values comparable to those reached in complete 
diagonalization studies [8|. 

Another exactly solvable class of quantum spin chains are spin- 1/2 XY chains. Rigorous 



analysis o 
fermions 



these systems is based on exploiting the Jordan- Wigner transformation to spinless 



Q. 



(For a relation between the Bethe ansatz method and the Jordan- Wigner 



approach for the spin-1/2 XXO (i.e. isotropic XY) chain see Ref. [iQj.) Although after 
applying the Jordan- Wigner transformation to the spin-1/2 XY chains one faces a system 
of noninteracting spinless fermions the calculation of the spin correlation functions is not a 
trivial problem because of the nonlocal character of the transformation. Thus, the zz spin 
correlations are related to the two-fermion (density-density) correlations, whereas, e.g., the 
xx spin correlations are related to many-fermion correlations. Accordingly, the zz dynamics 



is well studied 
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are rather scarce 
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^whereas closed-form expressions, e.g., for the xx dynamic quantities 
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151 ] (see also references in Ref. Il6l ). 



In the present paper we consider several quantum spin chains with regular alternation 
or random disorder in the nearest-neighbor interactions and follow the effect of such mod- 
ifications on the dynamic structure factors. The inhomogeneity introduced refers mainly 
to the sign of interactions and may mimic the ferromagnetic or antiferromagnetic types of 
nearest-neighbor exchange coupling. The interest in models of such a kind is not purely 
theoretical. Recently some organic and inorganic magnets have been recognized as alternat- 



ing sign [17(|, random bond [18| and alternating random bond [19|, |20|, |21] antiferromagnetic 
spin chains. The dynamic study of the quantum spin chain material with bond randomness 
BaCu 2 (Sii_ x Ge x ) 2 07, x = 0.5 using inelastic neutron scattering revealed that its dynamic 
structure factor can be fitted by the Miiller ansatz [l^ surprisingly well [3] . The correspon- 
dence between the dynamic properties of the random-bond Heisenberg antiferromagnetic 
spin chain and the BaCu2 (Sii-zGe^Or compound has been confirmed numerically by the 
quantum Monte-Carlo method [22j. 

In our calculation of dynamic quantities we use appropriate transformations to eliminate 
the inhomogeneity from the spin Hamiltonian arriving at the homogenous model the dynamic 
properties of which are well known. Thus we reduce the complex behavior of dynamic 
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quantities for periodic/random quantum spin chains to the known dynamic properties of 
the homogenous model. In what follows we deal with spin-1/2 isotropic XY (XX or XXO) 
chains since the dynamic quantities for the more general case of the XXZ Heisenberg 
exchange interaction are less known. 

The paper is organized as follows. To the end of this section we introduce the spin model, 
the quantities of interest and recall some results for the dynamic quantities obtained within 
the Jordan- Wigner fermionization approach which are used in the following sections. In 
Sec. [Til we consider the spin-1/2 XX chain with regularly alternating or random sign of 
the XX exchange interaction. In Sec. II III we consider the spin-1/2 XX chain with the 
Dzyaloshinskii-Moriya interaction the sign of which may either vary regularly along the 
chain (or it has a regularly varying component in addition to a constant component) or may 
acquire its sign randomly. We summarize our findings in Sec. IIV1 

We consider the following Hamiltonian of a one- dimensional spin s = 1/2 XX model 
with two-site interactions which can be examined rigorously within the framework of the 

CI 

Jordan- Wigner approach |9[: 

H = (j n (SnSn+l + Sn S n+l) + D n (s^Sn+l ~ S n S n+l) + ^ S n) ■ 

n 

Here J n is the exchange XX interaction between neighboring sites n and n + 1, D n is 
the ^-component of the Dzyaloshinskii-Moriya interaction between these sites, and Q is the 
external transverse (z) magnetic field. The sum in (II. ip runs over all N sites; the boundary 
conditions (periodic or open) are not essential for the quantities considered below which we 
calculate in the thermodynamic limit iV — > oo. 

We are interested in the dynamic structure factors of the spin model (11. ip (defined most 
conveniently for periodic boundary conditions, so that m = N is equivalent to m — 0) 

S^(K |W ) = -££exp(-iKm) / dtexpM) (^(t)^ +m ) - (^)(^ +m )) , (1.2) 

where a, (3 = x, y, z. These experimentally accessible quantities contain important in- 
formation about the spin model (11. ip . By symmetry arguments S xx (n, u) = S yy (K,u), 
S xy ( k, uj) — — Sy X ( — k^lo). Therefore, in what follows we may focus only on S xx (k, cj), 
S xy (n,u) and S zz (k,uj). Moreover, the model (II .ip implies that (s*) = (s y n ) = and hence 
the second term in the parentheses in Eq. (11.21) may be omitted if a, (3 = x, y. 
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Consider first a uniform chain (11.11) with J n = J, D n = 0. Again by symmetry argu- 
ments S zz (k, uj) is insensitive to a sign change of the exchange interaction J — » — J whereas 
S xx (k,u) — > S xx (k^tt,uj), S xy (K,ui) — > S X y(K^7T,uj). Next, from Refs. 

S zz (k,uj) = / dKxn Kl (1 - n K+Kl ) 5 (uj + A Kl - A K+Kl ' 
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we know that 
(1.3) 

where A K = Q + J cos k is the elementary excitation energy of the Jordan- Wigner fermions 
and n K = 1/(1 + exp (j3A K )) is the Fermi function. Obviously the zz dynamic structure 



factor (11.31) is governed by a continuum of two-fermion (particle-hole) excitations 12J. Let 
us introduce the following characteristic lines in the k-uj plane 
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where a = arccos (|f2|/| J|) varies from 7r/2 (when Q = 0) to (when \Q\ = |J|). The 
ground-state S zz (k,uj) is nonzero for \Q\ < \ J\ and in a restricted region in the k-uj plane 
(we assume \k\ < ir, uj > 0) with the lower boundary uji(k) = u^'(k) and the upper 



boundary uj u (k) 



uj 



W(k) if < U < 



7T 



2a or ^(ac) = u/ 3 )(/i;) if 



7T 



2a < \k\ < 



7i. Moreover, S zz (k,u) exhibits a finite jump (increasing its value by 2) along the middle 



boundary uj m (K) 



UJ 



(2), 



k), 7t — 2a < \k\ < 7r. Finally, S zz (k, u;) shows a van Hove singularity 



along the curve uj s (k) = uj( 3 \k). As temperature increases the lower boundary becomes 
smeared out and finally disappears. The upper boundary is given by uj^ 3 \k) and S zz (k,uj) 
becomes field- independent in the high-temperature limit. 

The xx /xy dynamic structure factor is governed by many-fermion excitations and there- 
fore is a m uch ^complicated ouantity (the two-fennion contribution to was 



discussed in Refs. |23|, |2J). However, the ground-state S xx (k } uj) and S xv (k,uj) can be easily 
calculated for strong fields \Vt\ > \J\ [lj]] 



7T 



S xx (k, uj) = isgn(£l)S xy (K } uj) = —5(uj— |0| — J cosk) . 



;i.s) 



Eq. (11. 5ft shows that all the spectral weight in this case is concentrated along the curve 



At sufficiently low temperatures (A^T/| J\ = 0.01 . . . 0.05) we know from numerics (see Ref. 



16l ) that although S xx (k,u>) and S xy (K,u) are not a priori restricted to a certain region in 



the k-oj plane (and indeed these quantities have nonzero values throughout the k-uj plane), 
nevertheless their values are rather small outside the two-fermion excitation continuum dis- 
cussed above. More precisely, the xx and xy dynamic structure factors show washed-out 
excitation branches roughly following the boundaries of the two-fermion excitation contin- 
uum (see Eq. ( 11.41) ) for J < or following these boundaries shifted along the K-axis by ir 
for J > 0. In the high-temperature limit we have k], [ljj] 



S xx (k,u) = -^rj [ exp [ - v ~ T2 "" 7 ) +exp 





iS xy {K,u) = ^ ( exp ( - v ~ T2 "" 7 ) - exp [ - v ~ ' T2 "" 7 ) ) , (1.7) 

i.e. the xx and xy dynamic structure factors in this case are ^-independent and display 
Gaussian ridges at u = ±Q. 

Similar results on the dynamic properties of the dimerized spin- 1/2 XX chain (i.e. with 
J n = J (1 — (— l) n <5), where < 5 < 1 is the dimerization parameter, and D n = in Eq. 



(11. ip ) can be found in Ref. |25| (and references therein). The dynamic properties of the 



uniform spin-1/2 XX chain with the Dzyaloshinskii-Moriya interaction (i.e. with J n = J, 
D n = D in Eq. (11.11) ) were discussed in Ref. [26. 

In what follows (Sec. [Til and Sec. IHIj) we use the results recalled here to examine the 
dynamic properties of quantum spin chains with special types of periodically varying or 
randomly distributed interspin interactions. 



II. SPIN-1/2 XX CHAIN WITH PERIODICITY/RANDOMNESS IN THE SIGN 
OF EXCHANGE INTERACTION 

In this section we consider the spin model with the Hamiltonian (11 .ip assuming J n = X n J 
with A n = ±1 and D n = 0, i.e. the exchange interaction between the sites n and n + 1 may 
be either antiferromagnetic if A n J > or ferromagnetic if X n J < depending on the given 
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sequence {Ai, . . . , \n}- Let us perform a gauge transformation 



S n ~> S n — ^1-^2 • • • A n _iS^, 
S l -»• S l = A 1 A 2 • • • K-lSl, 

S n ~^ S n — S n (2-1) 

after which the Hamiltonian H transforms into the Hamiltonian H of the homogeneous 
model with exchange constant J n = J (up to an inessential boundary term). (We denote 
the quantities related to the transformed (homogeneous) model by a tilde.) Obviously, 
according to (12.11) the zz dynamic structure (as well as all thermodynamic quantities) does 
not feel an inhomogeneous sequence of signs {Ai, . . . , Ajv}- In contrast, the xx and xy 
dynamic structure factors do depend on {Ai, . . . , Xn}- Below we consider separately the 
cases of periodic sequences and of random sequences of signs. 



A. Periodic case 



We begin with the case of period p = 2, i.e. {A n } = {1, —1, 1, —1, . . .}. After performing 
the transformation (12. ip we have §21-1 = s 2j-n $2j = ( — l)" ?+ls 2j; J = 1,2,... (here 

and to the end of the paper a, (3 = x, y) and therefore according to ( 11.2jl we can write 

1 N poo 

Sapfa, w) = ^'y] exp (-i«m) / dt exp (iut) a m (s"(t)s{ +m ) 
* 1 J -00 



m=l 

N 



1 . \ r°° 

-^exp(-iKm) / dtexp(iujt)b m (s%(t)4 +m ), (2.2) 
A m=i J -°° 



where {ai, a 2 , a 3 , . . .} = {1,-1,-1,1,1,-1,-1,...}, {61, b 2 , 63, • • •} = 
{-1,-1,1,1,-1,-1,1,...}. Noting that a m = ((1 - i)/2) exp(i7rm/2) + ((1 + 
i)/2) exp(3i7rm/2) and b m = a m+ i we immediately find from Eq. (I2.2p that 



(2.3) 



On the l.h.s in Eq. ( 12. 3p we have the dynamic structure factors for the periodic chain 
S a p(K,uj) whereas on the r.h.s. in Eq. (12.31) the dynamic structure factors S a p(K,uj) refer 
to the uniform chain with the exchange constant J; the latter quantities were discussed in 
Sec. [H These calculations can be easily extended for periodic chains of larger periods. For 
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FIG. 1: S xx (k,lu) for the spin-1/2 XX chain with periodic sequences of exchange interactions 
{J, — J, J, —J, . . .}, J = 1 (upper panel a) and {J, J, — J, J, J, — J, . . .} with J = 1 (middle panel b) 
and with J = — 1 (lower panel c) for f2 = 0.25 at low temperature, (3 = 20. 

example, for p = 3 with {A„} = {1, 1, — 1, 1, 1, —1, . . .} after performing similar calculations 
we arrive instead of Eq. (12. 3p at 



S a p(K, id) = ~£ a/ 3 (k + ^, LO ) + ^£ a/3 (/C + 7T, Ct») + ^ a/3 + W^j . 



(2.4) 



To illustrate the effect of a regularly alternating sign of exchange interaction on S xx (k, uj) 
we display this quantity calculated according to Eqs. (I2.3p . (12.41) in Fig. [H Evidently, in the 
high-temperature limit owing to ^-independence of S xx (k,uj) (see Eq. (II. 7p ) regular alter- 
nation of the exchange interaction signs does not manifest itself in the xx dynamic structure 
factor. However at low temperatures it may lead to rather intricate frequency /wave- vector 
patterns (see Fig. [p. Interestingly, we may reproduce the sequence {A n } knowing the num- 
ber of soft modes Kq and their position. In the limit T = and \Q\ > \J\ we may insert 

8 



Eq. (I1.5p into the r.h.s. of Eqs. (12. 3p . (12. 4p to find that the spectral weight is concentrated 
along the curves which follow from Eq. (jl.6p after corresponding shifts along the K-axis. 

The results of this subsection are complementary to the earlier results on thermodynamic 
and dynamic properties of periodic spin- 1/2 XX chains (see Refs. |27|, |28|, |29| and references 
therein) . 

B. Random case 

We now proceed with the case of randomly distributed signs of exchange interactions 
assuming {A n } to be a sequence of independent random variables each with the following 
bimodal probability distribution 

p(An) = P 5(X n + 1) + (1 - p)5(X n - 1), (2.5) 

where < p < 1. We are interested in random-averaged quantities and denote the average 



over all realizations of randomness as (...) = Y[ n dA n p(A n )(. . .). Random chains of that 
type (in fact, for more general XX Z coupling) were studied in Refs. 
Exploiting the gauge transformation (12. ip and Eq. (12.51) we find 



ad, 



31 



(snt>j + J = (1 - 2p) N <S?(*)^ + m>- (2-6) 

Introducing the correlation length £ = — l/ln|l — 2p\, the last expression (I2.6P can be 
rewritten as 

— — fexp(-M)( i?(i ) 3 1 + J, 0< P <i, 

( 1 t-ir«p(-f) i<p<i. ( 

As a result, the random- averaged dynamic structure factors (11.21) can be written as follows 
S a/3 (K,uj)= ^ exp (-mm - ^f] j dtexp (jut) (s°{t)s? +m ); (2.8) 

m=0,±l,±2,... ^ ^ ' 

here < p < 1/2. If 1/2 < p < 1, a factor (-l) m (see Eq. (EZD) should be taken into 



account in Eq. (12. 8p and the resulting expression S a p(K,u) for 1/2 < p < 1 corresponds to 



S Q /3( K T ^5 w ) i n formula (12.81) . We use Eq. (12. 8p to compute S a p(n^uo) through the known 
results for (s? (t)s f + m ) for the uniform chain with exchange constant J obtained analytically 



or numerically 



16 



32j (see Fig.[2]). 




-] -Cl5 ri Q_5 1 -1 -Cl5 ij QJ 1 

!. " t. 7i 



FIG. 2: S xx (k, uj) for the spin-1/2 XX chain with random-sign exchange interaction J = 1, 
SI = 0.01 (left panels, a ...e), fi = 0.25 (right panels, f ...j), (3 = 20. The values of p in Eq. 
(I23|) are as follows: from top to bottom p = 0.1 (f ss 4.48), p = 0.25 (£ ps 1.44), p = 0.5 (f = 0), 
p = 0.75 (f ~ 1.44), p = 0.9 (f « 4.48). 



Let us consider the case T = 0, \Q\ > \J\ when the xx dynamic structure factor S xx (k, uj) 
is given by Eq. (jl.5p . In the site representation we have 

dt exp (iut) (s^(t)s x j+m ) = — exp (i/cm) SU(«, u) 

i Z' 7 ' 

= — y dfi; exp (i/?m) 5 (a; — — J cos «) 

= T T {mK °\ o - + | J|) f + | J| - w) , (2.9) 
2 1 J sin Ko | 

where k = arccos((a;- |0|)/J). [In Eq. (Q and Eqs. (12101 . ([2~TT]) the ^-functions simply 
indicate the frequency range within which the equation co>— \Q\ — J cos k = has the solutions 
k = ±Ko-] Consider first the case < p < 1/2. After substitution of (12.91) into (12. 8p and 
some simple calculations one finds 



S xx (k,u) = ^ exp (-'mm - j S xx (m,u) 

m=0,±l,±2,... ^ ^ ' 

exp | — iftm — - — -J cos(m/«o) 

14-14-9 V ? / 



2 J sin Kn 

1 U| m=0,±l,±2,. 

x6(uj - \n\ + \J\)9(\n\ + \j\-uj) 

Y Jsinh j ^ J cosh | — (u — |Q|) cos 

2-y/ J 2 — (lu — |fi|) 2 ^u; — |fi| — J cosh I cosk^ + J 2 sinh 2 | sin 2 k 

x# (w - |n| + \J\)6 (|fi| + | J| - w) . (2.10) 



If 1/2 < p < 1, S xx (k, to) follows from Eq. (12.101) after the change k — > K=)=7r. One can easily 
note that Eq. (12.101) transforms into Eq. (11.51) in the nonrandom limit l/£ — > (i.e. p — > 
or p — > 1) [to show this one has to exploit the relation lim r ^ +0 (r/ ((a; — uj ) 2 + T 2 )) = 
ttS (uj — ujq)]. In the opposite limit of a completely random system l/£ — > oo (i.e. p — > 1/2) 
Eq. ( 12. 101) becomes 

1 



S xx (k,uj) = — =6(oj-\n\ + |J|)0(|fi| + \J\-w). (2.11) 

2yjJ 2 - (u- \n\) 2 

One immediately recognizes that Eq. (12. lip contains the contribution of only the autocorre- 
lation function (as it should be since the correlation leng th £ tends to zero) and since in the 
limit considered 4(s|(t)sJ) = (l/N) J^ K exp(— iA K t) 



151 ] the xx dynamic structure factor is 



proportional to the density of states of elementary excitations p(E) = (l/N) J2 K 5 (E — A K ) 



i.e. S xx (k,u) = (n/2)p(u) independent of k. For other values of p S xx (k,ou) (12. 10[) is 
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restricted to the frequency region \Q\ — \J\ < uo < + \J\ and shows square-root sin- 
gularities as to —>■ \Q\ ± |J|. The frequency profiles at fixed k resemble (although are not 
identical to) Lorentzian shapes centered at uj = \Q\ + J cosh (l/£) cosk with the line width 
T = | Jsinh(l/£) sin/c|. 

For nonzero temperature, T ^ 0, and for subcritical field values, \fl\ < J, (12.81) must be 
evaluated numerically (see Fig. [2]). In the case p = 1/2 the correlation length £ — > and 



one expects only the autocorrelation function to contribute to the ^-independent S xx (k,u>) 
and the frequency shape for any k is determined by the tu-dependence of S xx (0,u) = 
f™oo dt exp (iut) (Sj(t)5j) (/t-independent stripes near frequencies which dominate the au- 
tocorrelation function). 

We note some similarities to recent numerical results on the spin- 1/2 Ising chain in 

In particular, the horizontal (k- independent) stripe-like 



a random transverse field 



patterns in Fig. [2] resemble the results of Ref. |33| for strong disorder. This is to be expected 



since for strong enough disorder only local correlations survive and lead to a ^-independent 
dynamic structure factor. 

The scheme presented here can be also easily adapted to more complex models where 
alternation and randomness are mixed. For example, the ferromagnetic-antiferromagnetic 
random alternating quantum spin chain compound (CH3)2CHNH3Cu(Cl x Bri_ x )3 can be 
viewed as a spin-1/2 random alternating quantum Heisenberg chain [21( 

), (2-12) 

n 

where Jm-\ — J is the weak uniform exchange bond, Jm = ^2nJ is the strong random-sign 
exchange bond and {A2«} is the sequence of independent random variables each with the 
bimodal probability distribution ( 12. 51) . If we restrict ourselves to isotropic XY interactions 
between spins in ( 12.121) . the randomness can be excluded from the Hamiltonian by a slightly 
modified gauge transformation Sa n _ x = rim=i Sg n = s% n YY^i ^2m, n = 2,3... 
obtaining finally the Hamiltonian of a dimerized XX chain with the periodically varying 
exchange couplings J, 2 J, J, 2J. . .. The random-averaged dynamic structure factors can be 
calculated analogously to (I2.6P - (I2.10p . 
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III. SPIN-1/2 XX CHAIN WITH PERIODICITY/RANDOMNESS IN THE SIGN 
OF DZYALOSHINSKII-MORIYA INTERACTION 



We now consider the spin model with the Hamiltonian fll.il) assuming J n = J and D n = 
X n D with X n = ±1. (We note that the case J n = X n J, D n = D may be analyzed on the basis 
of the results reported below after exploiting the unitary transformation discussed in Ref. 



34].) It is generally known 



26 



35 



36l . 1371 . 1381 ] that the Dzyaloshinskii-Moriya interaction D r , 



can be eliminated from the Hamiltonian H (11. ip (up to an inessential boundary term) by 
the spin coordinate transformation 



COS ( 



s x n + sin <j) n s 



x 



y 
n > 



S n ~* S n = ~ sin 4>nS n + COS (j) n S 



V 

n ■ 



(3.1) 



where n = Y^m=o V 9 ™.' V'o is an arbitrary angle which is usually assumed to be zero and 
tamp m = D m /J, to = 1,2.... As a result, one faces the Hamiltonian H (11.11) without the 
Dzyaloshinskii-Moriya interaction, however, with a renormalized XX exchange interaction 
J n = sgn( J) *J J 2 + D\ . In the uniform case, when D n = D, the unitary transformation 
(13.11) was used in the recent studies of dynamics of quantum spin chains [26|, |37|, l38| . In 
this section we consider separately the two cases of periodically varying Dzyaloshinskii- 
Moriya interaction and of random-sign Dzyaloshinskii-Moriya interaction focusing on the 
xx dynamic structure factor S xx (k, uj). 



A. Periodic case 

We begin with the case p = 2 with {A„} = {1, -1, 1, -1, . . .}, i.e. D n = (-l) n+1 D. Then 
we have to put in Eq. (13.11) <p m = (— l) m+1 <p, ip = arctan(D/ J). Moreover, it is convenient 
to assume tpo = —f/2. Then <p n = (—l) n ip/2 and the inverse transformation to the one 
given by (13. 1 1) reads 

< = cos|^-(-l) n sin|^, 
^=(-l)"sin|5: + cos| S X 

S n = *v (3-2) 
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By substituting Eq. (I3.2p into Eq. (11.21) one immediately finds that the xx dynamic structure 
factor S xx (k, uj) of the XX chain with the alternating Dzyaloshinskii-Moriya interaction 
D, —D,D, —D, . . . can be expressed through the xx dynamic structure factor S xx (k, uj) of 
the uniform chain with only XX exchange interaction J = sgn( J)\/ J 2 + D 2 as follows 

S xx (k, uj) = cos 2 ^ S xx (k, uj) + sin 2 ^ S xx (k ^ tt, uj). (3.3) 

We notice here that in the case when in Eq. (11.11) J n = J, D n = Dq, that is, for uni- 
form XX and Dzyaloshinskii-Moriya couplings, the relation for S xx (k, uj) is quite different: 
S xx (k,uj) = (S xx (K-(p,uj) + S xx (K + (p,u) + iS xy (K-(p,u)-iS xy (K + ip,u))/2] here S a/3 (K,uj) 
is related to the uniform chain with only XX exchange interaction J = sgn(J)^y J 2 + D 2 , 



(see Ref. |26| ). It is worth therefore to consider also the more complicated case of the chain 
(II. ip with J n = J and D n = D Q — (—l) n D. This choice of a dimerized Dzyaloshinskii-Moriya 
interaction covers both limiting cases (i) of the alternating-sign Dzyaloshinskii-Moriya inter- 
action when Dq = and (ii) of the constant Dzyaloshinskii-Moriya interaction when D = 0. 
Exploiting the transformation (13. ip with (p m = arctan((D — (—l) m D) /J) we arrive at a 
chain without the Dzyaloshinskii-Moriya interaction but only with the dimerized XX ex- 
change interaction J n = sgn(J)y^ J 2 + (Dq — (—l) n D) 2 . To find the relation between the 
xx dynamic structure factor S xx (k, uj) of the XX chain with the dimerized Dzyaloshinskii- 
Moriya interaction and the dynamic structure factors S a f3(n, uj) of the dimerized XX chain 
without the Dzyaloshinskii-Moriya interaction we proceed as follows. First, we note that 
exploiting (ETTJ) in Eq. flL2D yields 

S xx (k, uj) = — exp (— mm) / dt exp (iut) 

N j=l m=l 

x (cos {<f> j+m - 0,) (s*(t)s x +m ) - sin {<jj j+m - 0,) (s x (t)s y j+ J) . (3.4) 

After introducing the notations <£> = arctan((.D + D)/J), ip c = arctan((D ~ D)/ J) and 
ip = ((f ± <f c )/2 we can write (f)j +m — <f)j = rrup + + (— 1) J (((— l) m — l)/2) tp~. Then after 
inserting this result into Eq. (13.41) and some manipulations Eq. (13.41) becomes 

S xx (k, uj) = - cos 2 ^- (J3 xx {k - (p + , uj) + S xx (k + y? + , uj) 
+iS xy (K - Lf + , uj) - iS xy (n + tp + , uj) 

+ \ Sin2 {^ XX ( K ^ 71 ~ tf+J ^ + ^ XX ( K ^ n + ^ 

+iS xy (K =F tt - v? + , uj) - \S xy {K =F vr + ip + , uj)j . (3.5) 
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FIG. 3: S xx (k, u) for the chain (jl.ip with J n = 1, D n = Dq — (—l) n D, Dq = (upper panel a), 
Dq = 0.25 (lower panel b), D = 0.5, ft = 0.25 at low temperature, = 20. 



Eq. (13.51) in the limit Dq = transforms into (13.3j) since ip + = 0, 9? = 9? = arctan(D/ J). 
Eq. (13.51) also contains the result of Ref. |26| in the limit D = since <p + = arctan(Do/J), 
= 0. 

In Fig. [3] we illustrate the effect of the dimerized Dzyaloshinskii-Moriya interaction on 
the xx dynamic structure factor at low temperatures. The panel a corresponds to the case 
D = (S xx (k,uj) is obtained using Eq. (13. 3p ). whereas the panel b corresponds to the case 
Dq 7^ (S xx (k, uj) is obtained using the more general Eq. ( 13. 5j) ) . 



B. Random case 



Finally, we pass to the case when the Dzyaloshinskii-Moriya interaction D n = X n D is 
given by a sequence of independent random variables {X n } each with the bimodal proba- 
bility distribution H 2 . 5 1) . For a specific realization of the signs of the Dzyaloshinskii-Moriya 
interaction we can eliminate D n from the Hamiltonian H (11.11) by the transformation (13. ip 
with ip m = \ m (p, if = arctan(/}/ J) arriving at the model H with only XX exchange interac- 
tion J = sgn( J)\/ J 2 + D 2 . To calculate the random- averaged xx dynamic structure factor 
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we need 



<*f (*K+J = cos {4> j+m - 0i)(s|(t)SJ +m ) - sin (<P 3+m - 0,)(SJ(t)^ +m ). (3.6) 
Noting that 

cos((Ai + • • • + A m ) <p) = - {pexp (-iif) + (1 -p)exp {if)) m 

+^ (pexp (itp) + (1 -p) exp (-i<^)) m 

m 

= (cos 2 <£> + (1 — 2p) 2 sin 2 2 cos (m arctan ((1 — 2p) tany?)) , 



sin ((Ai + . . . + \ m ) if) = (cos 2 cp + (1 — 2p) 2 sin 2 y?) 2 sin (m arctan ((1 — 2p) tan<£>)) (3.7) 

and introducing the notations £d = — 1/ In a/ cos 2 <p + (1 — 2p) 2 sin 2 <p, <pu = 
arctan ((1 — 2p) tamp) one finds that 

<sj(t)sj+j = exp f-^J (cos (m^) (sj(t)s* +m ) - sin (m<^) (SJ(t)*?+m» ■ ( 3 - 8 ) 
Using Eq. (I3.8P the random-averaged xx dynamic structure factor can be written as follows 
S xx (k, uj) = exp ( — inm — J f dtexp(iut) 

m=0,±l,±2,... ^ ^ D ' J' 00 

x (cos (m<p D ) (s*(t)s* +m ) - sin (m<p D ) (§ x 3 (t)s y j+ J) . (3.9) 
On the r.h.s. in Eq. (13.91) we have the correlation functions of the uniform XX chain with 



the exchange constant sgn( J) \J J 2 + D 2 . We use Eq. (13.91) to calculate S xx (k, uj) for the 
model with the random-sign Dzyaloshinskii-Moriya interaction through the known results 
for (s"5j +m ) la, [32]. The results are shown in Fig. HI The correlation length £p attains 
its minimal (nonzero) value — 1 / In | cos ip\ at p = 1/2. A comparison of Figs. Hb andHh to 
Figs. Ek, EF shows that S xx (k, uj) looks very similar for weak disorder (p = 0.1) in the sign 
of the exchange interaction on one hand and for maximum disorder (p = 1/2) in the sign 
of the Dzyaloshinskii-Moriya interaction on the other hand. It looks as if the random-sign 



Dzyaloshinskii-Moriya interaction at p = 1/2 does not manifest itself in S xx (k, uj). This 
similarity becomes evident if we notice that the correlation lengths £ and £d are of the same 
order for the considered conditions (see the captions to Figs. El HJ) and tpo tends to zero 
which cancels any signals of the Dzyaloshinskii-Moriya interaction due to formula (13.91) . 



We can proceed with analytical calculations for the case T = 0, > \J J 2 + D 2 . Com- 
paring Eq. (I3.9P and Eq. (I2.10p we see that the result we are interested in follows from 
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FIG. 4: S xx (k, uj) for the chain (11.11) with J n = 1, D n = X n D, where A n is a random variable with 
the probability distribution (12.51) . D = 0.5, = 0.01 (left panels, a . . . e), O = 0.25 (right panels, f 
. . . j), /3 = 20. The values of p in Eq. (|2.5p are as follows: from top to bottom p = 0.1 (£d 26.77), 
p = 0.25 ~ 12.31), p = 0.5 « 8.96), p = 0.75 (£ D » 12.31), p = 0.9 » 26.77). 



ar, 



Eq. (12.101) after the changes J — * J = sgn( J)\J J 2 + -D 2 , k — > k — sgn(fi)<^£>. In particu 
for the nonrandom case when p = or p = 1 we recover the result reported in Ref. 
S xx (k,uj) = (ir/2)S — |0| — Jcos(k — sgn(f2)<^)J. 

Interestingly, we can extend the scheme explained above to more complicated random 
chains assuming J n = J cos (/(A n )), D n = J7"sin (f(X n )), where f(x) is an arbitrary function, 
for example, f(x) = A + Bx, and A n is a random variable with an arbitrary probability 
distribution p(X n ) (not necessarily with the bimodal probability distribution (12.51) ). After 
exploiting the transformation (13.11) with ip m = A + BX m we arrive at the Hamiltonian H 



given by Eq. (II. ip with J n = J D n = 0. (sj(t)sj +m ) is again given by Eq. (13.61) . however, 
Eq. (13. 7p now reads 



cos (mA + B (Ai + . . . + A m )) = \F(B) \ m cos (m (A + arg(F(5)))) 



sin(mA + 5(A 1 + ... + A m )) = \F(B)\ m sin (m (A + axg(F(B)))) , (3.10) 

where 

F(5) = y dA n p(A„)exp(i J BA„) = \F(B)\ exp (i arg(F(S))) (3.11) 

is the characteristic function of the random variable A n . Now we introduce the notations 
£d = -l/ln|F(S)|, <^d = A + arg(F(5)) and arrive at Eq. (EEJ) and Eq. (ET5]I . For the 
model with the bimodal distribution considered earlier we have to put J = sgn( J)y/ J 2 + D 2 , 
A = 0, B = arctan(D/J) = <p, and therefore F(B) = pexp(-iip) + (1 — p) exp (iyj), 



|F(_B)| = y cos 2 + (1 — 2p) 2 sin 2 ip, axg(F(B)) = arctan ((1 — 2p) tan<^) and we reproduce 
Eqs. (I3.8p . (13.91) with the expressions for ^ and tpo given just before Eq. (13.81) . 



IV. CONCLUSIONS 



To summarize, we have considered a number of inhomogeneous (periodic or random) 
spin-1/2 XX chains, to examine their dynamic properties. The models considered are 
distinguished by the possibility to eliminate the inhomogeneity from the spin Hamiltonian 
by a suitable unitary transformation (see Eqs. (12.11) . (I3.ip ) and therefore to reduce the 
problem to the well known one for the uniform model. We use exact analytical and precise 
numerical data to analyze the dynamic structure factors of the periodic/random spin-1/2 
XX chains. The models considered show rather complex behavior which, however, can be 
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explained by the corresponding properties of the basic uniform model. Thus, for the periodic 
chains only the correspondingly modified characteristic curves Eqs. (11.41) . (jl.6p are seen in 
the complex pattern displayed by the dynamic structure factor at low temperatures. In 
the high-temperature limit only Eq. (11.71) is relevant. In the cases considered the observed 
complexity has a simple origin. We also stress here that we have reported rigorous analytical 
results for dynamic structure factors of some periodic/random quantum spin chains. In 
comparison, direct numerical treatment of random quantum spin chains would imply many 
calculations of dynamic quantities for different realizations of the random couplings and 
a subsequent average over these realizations, which altogether would require an enormous 
amount of computer time. 

It is interesting to note that the effects of temperature and of random couplings on 
the xxjxy dynamic structure factors are different (compare Eq. (II. 7p and Eqs. ( 12.81) . 
(12. lip ). Although in both cases only the autocorrelation function determines the dynamic 
structure factor (for sufficiently high temperature or sufficiently strong randomness), at high 
temperatures the dynamic structure factor is ^-independent and shows Gaussian rid ges (see 



Eq. fll.7p ). That is due to the Gaussian time decay of the autocorrelation function [13l. Il4] 
which should be contrasted to the slow long-time decay of the autocorrelation function at 
low temperatures. 

The spin chain models discussed in our study are obviously of a rather special kind, 
and it would be highly desirable to obtain reliable results also for more general types of 
inhomogeneity in the interspin couplings, where not only the signs but also the absolute 
values of the couplings vary. For those more general models, however, the present methods 
are not applicable, and different methods or approximations have to be employed, such as 



in Refs. 
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for example. The special models treated in our present study will then be 
useful in providing a testing ground for the more general (but possibly less reliable) methods 
capable of dealing with a broader class of systems. 

Finally, the dynamic structure factors provide benchmarks for determining interspin in- 
teractions. In our paper we have demonstrated by some examples how periodic modulations 
or random variations in the signs of nearest-neighbor interactions manifest themselves in the 
dynamic structure factor. We note that the techniques used here for XX chains may also 
be applied to study dynamic structure factors of more general XX Z chains with periodic or 
random sign changes in the XY part of the interactions, provided that sufficiently precise 
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data for the corresponding uniform systems become available. 
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